Computational Study of Performance Characteristics for Truncated Conical 


Aerospike Nozzles 


Prasanth P Nair’, Abhilash Suryan’ and Heuy Dong Kim? 


1. Department of Mechanical Engineering, College of Engineering Trivandrum, Trivandrum, 695016, India 


2. Department of Mechanical Engineering, Andong National University, Andong, 760749, Korea 


Aerospike nozzles are advanced rocket nozzles that can maintain its aerodynamic efficiency over a wide range of 


altitudes. It belongs to class of altitude compensating nozzles. A vehicle with an aerospike nozzle uses less fuel at 


low altitudes due to its altitude adaptability, where most missions have the greatest need for thrust. Aerospike 


nozzles are better suited to Single Stage to Orbit (SSTO) missions compared to conventional nozzles. In the cur- 


rent study, the flow through 20% and 40% aerospike nozzle is analyzed in detail using computational fluid dy- 


namics technique. Steady state analysis with implicit formulation is carried out. Reynolds averaged Navier-Stokes 


equations are solved with the Spalart-Allmaras turbulence model. The results are compared with experimental 


results from previous work. The transition from open wake to closed wake happens in lower Nozzle Pressure Ra- 


tio for 20% as compared to 40% aerospike nozzle. 


Keywords: aerospike nozzle, plug nozzle, advanced rocket nozzles, thrust coefficient, base pressure. 


Introduction 


Conventional bell nozzles are not suited for Single 
Stage To Orbit (SSTO) missions on account of the fixed 
geometry which cannot compensate performance at va- 
rying altitudes. The performance is reduced at low alti- 
tudes due to overexpansion and at high altitudes due to 
under expansion. At low altitude where the atmospheric 
pressure is high, the flow gets compressed radially in- 
wards towards the nozzle exit. This can cause formation 
of Mach disc and flow separates from the boundary wall 
within the nozzle and hence, decreasing the efficiency as 
shown in Fig. 1(a). At higher altitudes, the atmospheric 
pressure decreases, causing the exhaust gases to continue 
to expand beyond the nozzle exit. This external expan- 
sion does not exert force on the nozzle wall, resulting in 
reduced thrust and hence, again decreasing the engine 
efficiency as shown in Fig. 1(b). There are studies on 
different types of advanced unconventional rocket noz- 
zles which can be used efficiently in the future to replace 
conventional bell nozzles. The different types of rocket 
nozzles discussed by Hagemann [1] are plug nozzle, ex- 
pansion-deflection nozzle, dual-bell nozzle, and 


dual-expander nozzle. All these nozzle concepts have 
altitude compensating capabilities. 

Plug (spike) nozzle can be described as bell nozzle inside 
out. Unlike in bell nozzle, the exiting gases in plug nozzles 
are not contained by the outer boundary. The spike prevents 
turbulent mixing in the exhaust flow that would occur in its 
absence, while maintaining isentropic expansion. For aero- 
spike nozzles, the spike is truncated which keeps the flow in 
closed wake rather than in open wake at high altitudes. 

At low altitudes, when the atmospheric pressure is 
high, the exhaust remains close to the nozzle contour, 
and at higher altitudes, recompression shock waves 
occur inside the plume, forcing the plume to remain 
column-shaped in a closed-wake state, irrespective of 
expanding exhaust plume shape. This maintains effi- 
ciency because the exhaust gases are not expanding 
past the physical nozzle, but still exerting force on the 
nozzle as shown in Fig. 2. The performance loss due to 
truncation at high altitudes is countered by providing 
secondary flow into the base region, hence the name 
aerospike as given by Rocketdyne. Aerospike or trun- 
cated plug nozzle was first proposed by Griffith [2] in 
1954 and then upheld by Rocketdyne from the beginning 


of 1960s [3]. Different methods for contouring of nozzle 
surface were proposed by Rao [4] and Angelino [5] to 
optimize the thrust. 
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Fig. 1 Flow phenomena for a conventional rocket nozzle: 


(a) Overexpansion; (b) Under expansion [1] 


An analytical method was presented by Migdal [6] for 
the design of axisymmetric annular nozzles to yield uni- 
form and axial flow. Sule et al. [7] studied base pressure 
characteristics of annular truncated plug nozzle. 

In 1997, Rocketdyne, conducted the Linear Aerospike 
SR-71 Experiment (LASRE) program to demonstrate the 
performance of the linear aerospike on a lifting body to 
simulate the X-33. This was accomplished by mounting a 
20% scale, X-33 fore body onto an SR-71 [8]. Tomita et 
al. [9-11] performed experimental studies on linear and 
clustered plug nozzles. Research was also conducted in 
Europe on clustered plug nozzles [12]. Bui et al. [13] 
conducted flight research on aerospike nozzle using solid 
motor rocket. 

In the current study, computational fluid dynamics 
(CFD) simulation is performed on a 15° annular conical 
aerospike nozzle (design Mach number 2.0) without free 
stream flow [14]. The simulation is conducted on 40% of 
the spike length and 40% of the spike length with differ- 
ent Nozzle Pressure Ratios (NPR). NPR is the ratio of jet 
stagnation pressure (P,;) and ambient pressure (P,). The 
numerical simulation results and experimental results are 
compared. 
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Fig. 2 Flow phenomena of an aerospike nozzle (a) Overex- 
pansion (b) Under expansion 


Experimental details and modeling 


Experimental study of Verma [14] is considered for 
validating the numerical results. Experiments were con- 
ducted in a 0.5m base flow facility, in a special purpose 
blow down-type tunnel. The results in the experiment 
indicated good mean flow uniformity and axisymmetry 
of the jet flow in the tunnel. 

Design Mach number 2.0 conical (15° half-angle) an- 
nular aerospike model was mounted on the central cylin- 
drical inner body. The truncated nozzle had two static 
pressure ports in the base to acquire base pressure. Exit 
nozzle radius (r.) was fixed as 25 mm. The annular gap at 
throat (h) is 9 mm, as shown in Fig. 3. Jet stagnation 
pressure P,; was changed to vary NPR [14]. 
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Fig. 3 Schematic of the experimental annular conical aero- 
spike nozzle model 


Numerical method 


Numerical modeling and analysis of flow for aero- 
spike nozzle with different plug shapes was done using 
ANSYS Fluent. Reynolds averaged Navier Stokes 
(RANS) solver was used along with one equation Spa- 
lart-Allmaras turbulence model. 

In this study, 2-D computation is done assuming 
axisymmetric flow. The domain varies from 60D (where 
D is exit nozzle diameter) horizontally to 10D vertically 
and 5D in the upstream of the throat as shown in Fig. 4 
(a). Inlet is at 2.5D from the upstream of throat. The 
mesh is generated in ICEM CFD, with the first grid point 


Table 1 Grid independence study for nozzle with 40% spike 


SL No. Cells Mass flow rate at nozzle exit 


Velocity at nozzle exit 


from either the spike or the cowl wall giving a y+.The 
mesh is finer at the throat region as shown in Fig. 4 (b). 

The grid independence study was conducted on the 
40% aerospike nozzle to reduce influence of grid on the 
computational result. The coefficient of thrust, mass flow 
rate, velocity and pressure at the exit were compared to 
select the computational grid as shown in Table 1. On 
comparison, the values are almost constant for coarse and 
fine mesh. Hence coarse mesh was selected for the 
computational study to reduce the computational time. 
The computational grid selected was 0.156893 million 
for 40% aerospike nozzle and 0.155133 million for 20% 
aerospike nozzle. The governing equations used are as 
follows: 
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Fig. 4 Domain of nozzle numerical model: (a) overall com- 
putational domain (b) closer view of computational 
3rid in the area of the nozzle region 


Pressure at nozzle exit Coefficient of Thrust 


(kg/s) (m/s) (Pa) 
1. 156893 0.682836 433.3473 67553.96 0.810273 
2. 310613 0.682829 433.3339 67558.63 0.810349 
3. 614188 0.682773 433.3197 67559.49 0.810236 


Mass Conservation Equation: 
For 2D axisymmetric geometries, the continuity 
equation is given by: 
Ps (pv,)+Z(pu,)+e=0 A) 
where x is the axial coordinate, r is the radial coordinate, 
vx is the axial velocity, and v, is the radial velocity. 
Momentum Conservation Equation: 


For 2D axisymmetric geometries, the axial and radial 


momentum conservation equations are given by: 
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Energy Equation: 

Turbulent heat transport is modeled using the concept 
of Reynolds' analogy to turbulent momentum transfer. 
The energy equation is thus given by the following: 


ô ô 
—(pE)+—|u;(pE+ p)] 
ôt Ox; (4) 


Where kr, is the thermal conductivity, E is the total 


energy, and (tj P is the deviatoric stress tensor. 


Spalart-Allmaras Turbulence Model: 

The transported variable in the Spalart-Allmaras 
turbulence model [15], © is identical to the turbulent 
kinematic viscosity. Transport Equations for Dis given 
by: 
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In these equations, G 


» is the production of turbulent 


viscosity and Y, is the destruction of turbulent viscosity. 


Og and C} are constants and v is the molecular 


kinematic viscosity. 


Comparison of flow structure with experimental 
result 


For 40% truncated spike, the Schlieren images of ex- 
perimental results of NPR 2.57, 4.15 and 7.0 are com- 
pared with shadowgraph obtained from numerical results. 
Due to truncation, trailing shock and recirculation occurs 
at the base region of nozzle. Flow features that can be 
observed from these results are compression shocks, ex- 
ternal jet boundary which developed from the cowl exit, 
the internal nozzle overexpansion shock originating from 
cowl lip, and the overexpansion shock formed on the 
nozzle surface. At low NPR, interaction between over- 
expansion shocks on nozzle surface results in flow sepa- 
ration. While increasing the NPR, the angle of these 
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shocks is changed as the internal nozzle operates from 
over expansion condition to under expansion condition. 
The flow conditions on the nozzle can be classified into 
three types as ‘A’, ‘B’ and ‘C’. In type ‘A’ flow, the 
change of angle for the overexpansion shocks from the 
internal nozzle alters the interaction between overexpan- 
sion shocks at the cowl lip and nozzle surface in a way 
that the reflected shock causes flow to reattach on the 
plug surface causing a separation bubble as shown in 
Fig. 5, which increases the pressure at the base. The reat- 
tachment of flow can cause undesirable side force. The 
flow on the plug further expands around the corner, 
which further decreases the base pressure. In type ‘B’ 
flow, the overexpansion shock stays exactly at the corner 
of the truncation as shown in Fig. 6. 

When internal nozzle begins to operate in the under- 
expanded condition, the expansion fan from the cowl lip 
inflicts the inner shear layer, this leads to the generation 
of a strong trailing shock that is visible in Fig. 7 for NPR 
7.0, this flow is type ‘C’, the nozzle operates in closed 
wake condition [14]. 
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Fig. 5 Flow structure for 40% spike at NPR 2.57: (a) Experi- 
mental Schlieren image [14]; (b) Numerical shadow- 
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Fig. 6 Flow structure for 40% spike at NPR 4.15: (a) Experimental Schlieren image [14], (b) Numerical shadowgraph 
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Fig. 7 Flow structure for 40% spike at NPRe 7.0: (a) Experimental Schlieren image [14], (b) Numerical shadowgraph 


Pressure ratio along the base for 40% spike 


The base pressure comparison of experimental and 
numerical results in 40% spike nozzle is plotted for va- 
rying NPR. The ratio of pressure at base (P,) and jet 
stagnation pressure (P,;) is plotted against NPR in Fig. 8. 
Two major flow regimes, open wake and closed wake 
can be identified [10, 16-18] from Fig. 8. The one in 
which base pressure ratio (P,/P,j) continues to decrease 
sharply and later, in which the base pressure ratio is con- 
stant. The open wake is based on the influence of am- 
bient pressure on base pressure. While analyzing Fig. 8, 
there is a decrease in base pressure ratio P,/P,j, till NPR 
5.6 for experimental result and NPR 6 for numerical re- 
sult, which corresponds to open wake regime, after the 
transition from open wake to closed wake the base pres- 
sure remains constant since base pressure is not affected 
by ambient pressure [16]. 


Performance characteristics: 40% spike 


Coefficient of thrust obtained in experimental study is 
compared with current numerical results with respect to 
varying NPR. The comparison for the nozzle with 40% 
spike is shown in Fig. 9. 
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Fig. 8 Comparison of experimental and numerical base pres- 
sure measurements for 40% spike; Py/Poj vs NPR 


——1m— Experimental 
—_4— Numerical 


o 
e 
Q 
0 2 4 6 8 10 12 
NPR 
Fig. 9 Comparison of experimental and numerical coefficient 


g = of thrust for 40% spike, C(x) vs NPR 


The thrust of the nozzles were found out using thefol- 
lowing formula [14, 19]: 
F =([mV, + (Parit — Pa ) Aerie | 


exit 
+ [(P, - P, )cos75dA + (P, — P, )Ay 


Where dA, is the elemental area on the spike surface, Pw is 
the wall pressure and P,, the ambient pressure. The 
coefficient of thrust of a nozzle is given by [19]: 

C,(x)=F/(4,Py) (7) 

The numerical Cris in match with the experimental 
result for 40% spike. Numerical results of 40% spike 
nozzle are used to compare the base pressure and C; of 
20% spike nozzle. For 20% spike nozzle, the base pres- 
sure decreases and becomes constant at NPR 5.7 as 
shown in Fig. 10. On comparing 20% spike nozzle with 
40% spike nozzle it is seen that transition from open 
wake to close wake regime happens at lower NPR. 
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Fig. 10 Comparison of base pressure measurement between 
20% aerospike nozzle and 40% aerospike nozzle, 
P,/P,j vs NPR 
Although truncation of nozzle leads to reduction in 
thrust as shown in Fig. 11, but for practical purpose 
transition at low NPR is preferable and at higher NPR 
larger base area is advantageous as it contributes in 
increased thrust [20]. The thrust loss can be compensated 
by providing a secondary flow through the base. 
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Fig. 11 Coefficient of thrust; comparison between 20% spike 
and 40% spike; C(x) vs NPR 


Conclusion 


The base pressure, shock pattern, flow structure and 
performance characteristics of numerical result agrees 
with the experimental data. The flow separation due to 
restricted shock separation at low NPR can result in 
undesirable side or lateral forces. In order to reduce this 
effect, truncated plug nozzles are suggested in which 
there is no reattachment of flow due to shorter length. 
Due to truncation there is the reduction in thrust but 


transition from open wake to closed wake occurs in 
lower NPR. 
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